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Recently, magnons, which are quasiparticles describing the collective motion of spins, were found 
to undergo Bose-Einstein condensation (BEC) at room temperature in films of Yttrium Iron Garnet 
(YIG). Unlike other quasiparticle BEC systems, this system has a spectrum with two degenerate 
minima, which makes it possible for the system to have two condensates in momentum space. 
Recent Brillouin Light scattering studies for a microwave-pumped YIG film of thickness d = 5 ^m 
and field H = 1 kOe find a low-contrast interference pattern at the characteristic wavevector Q of the 
magnon energy minimum. In this report, we show that this modulation pattern can be quantitatively 
explained as due to non-symmetric but coherent Bose-Einstein condensation of magnons into the 
two energy minima. Our theory predicts a transition from a high-contrast symmetric phase to a 
low-contrast non-symmetric phase on varying the d and H, and a new type of collective oscillations. 

PACS numbers; 75.10.-b, 75.60, 75.70, 75.85 



Bose-Einstein condensation (BEC), one of the 
most intriguing macroscopic quantum phenomena, 
has been observed in equilibrium systems of Bose 
atoms, like ''He [HE], ^^Rb [3] and ^^Na j4j. Recent 
experiments have extended the concept of BEC to 
non-equilibrium systems consisting of photons [5] 
and of quasiparticles, such as excitons [Bl , polari- 
tons [THl] and magnons [TUl E]. Among these, 
BEC of magnons in films of Yttrium Iron Garnet 
(YIG), discovered by the group of Deniokritov [11 - 
[T7] . is distinguished from other quasiparticle BEC 
systems by its room temperature transition and 
two-dimensional anisotropic properties. In par- 
ticular, the spin-wave energy spectrum of a YIG 
film shows two energetically degenerate minima. 
Therefore it is possible that the system may have 
two condensates in momentum space [18] . An ex- 
periment by Nowik-Boltyk et al. |17| indeed shows 
a low-contrast spatial modulation pattern, indicat- 
ing that there is interference between the two con- 
densates. Current theories [I9H24] do not describe 
the appearance of coherence or the distribution of 
the two condensates. 

This report points out that a complete descrip- 
tion of BEC in microwave-pumped YIG films must 
account for the 4th order interactions, includ- 
ing previously neglected magnon-non-conserving 
terms originating in the dipolar interactions. The 
theory explains not only the appearance of coher- 
ence but also quantitatively explains the low con- 



trast of the experimentally observed interference 
pattern. Moreover, the theory predicts that, on 
increasing the film thickness d from a small value, 
there is a transition from a high-contrast symmet- 
ric phase S for d < dc, with equal numbers of con- 
densed magnons filling the two minimum states, to 
a low-contrast coherent non-symmetric phase NS 
for d > dc, with different numbers of condensed 
magnons filling the two minimum states. In com- 
paratively thin films (d < 0.2/im) the same transi- 
tion can be driven by an external magnetic field H. 
For d > d* , where d* is another critical thickness 
{d* > dc), the sum of phases of the two conden- 
sates changes from tt to 0; for d = d* the system 
is in a completely non-symmetric phase with only 
one condensate, for which there is no interference. 
In the experiment of Ref.[T7j the film thickness d 
exceeded d* . We suggest that the phase transitions 
may be identified by measuring the contrast of the 
spatial interference pattern for various d and H. 
We also predict a new type of collective magnetic 
oscillation in this system and discuss the possibil- 
ity of domain walls in non-symmetric phases. 

Results 

Phase Diagram. — We consider a YIG film of 
thickness d with in-plane magnetic field H (see in- 
set of Fig. [T]). The 4-th order interaction of con- 
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densate amplitudes reads 

+ C[c'qCqCqC-Q + c'^qC-qC^qCq + h.c.]. (1) 

Here c±q and Cj_g are the annihilation and cre- 
ation operators for magnons in the two conden- 
sates located at the two energy minima (0, ±Q) in 
the 2-D momentum space (see. Fig[l]). The cocffi- 
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FIG. 1. The magnon spectrum in the kz direction for 
d = 5 /^m and H — 1 kOe. The inset is a schematic 
diagram of YIG film. 
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with ai ^ u"^ + Au^v^ + a2 2v?v'^ and — 
3uv{u'^ + w^). Here, u and v are the coefficients of 
Bogoliubov transformation (see the Methods sec- 
tion for details). S = 14.3 is the effective spin, N is 
the total number of spins in the film, M is the mag- 
netization and huj]\i = jAnM with gyromagnetic 



to the exchange constant and Fk — (1 



)/kd. 



Similar results for the coefficients A and B were 
obtained in Ref. jl9j. The coefficient C, which vi- 
olates magnon number conservation, has not been 
considered previously. Below we show that C is the 
only source of coherence between the two conden- 
sates. The three coefficients A, B and C, whose 
values as functions of H are shown in Figj2] for 
two typical values of d, determine the distribu- 
tion of condensed magnons in the two degenerate 
minima. Ref. |19| assumed a symmetric phase with 
condensed magnons in both minima having equal 
amplitudes and equal phases. Later, Ref. [20] as- 
sumed filling of only one minimum. More recently 
Ref. [23] considered Josephson-like oscillations by 
starting from two condensates with equal numbers 
of magnons but different phases. Our theory pre- 
dicts coherent condensates and the ratio of their 
amplitudes with no additional assumptions. 

In terms of condensate numbers N±q and 
phases the condensate amplitudes are c±q — 

N±Qe'^'^^ . Substituting them into eq.|l]) we find: 

V4 ^A{N% + N'i 



+2Ccos$(iv|iV5^ + nI^NIq). (3) 

Clearly the dipole energy depends on the total 
phase $ = 4>+ + </>-• To minimize this energy, 
for C > we have $ = tt and for C < we have 
$ = 0. Fig[2] shows that the sign of C changes for 
different d and H , which indicates a transition of 
$ between and tt. For both C > and C < 
the dipole energy establishes a coherence between 
the two condensate amplitudes. In contrast to a 
Josephson-like interaction, the sum rather than the 
difference of the two condensate phases is fixed. 

Since the total number of condensed magnons 
Nc — Nq + -/V_Q is uniquely determined by the 
pumping (see Methods) , the energy has only a sin- 
gle free variable, the so far unspecified difference 
5 — Nq — N^Q. In terms of Nc and 5 the conden- 
sate energy eq.([3]) is: 

^ \{A + B)Nl -{B- A)5'^ 



-2\C\Nc^N^-S^ 



(4) 



The ground state of the condensates depends on 
the criterion parameter A, defined as 

A = A~B+\C\. (5) 
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FIG. 2. Interaction coefRcients A, B and C (in units of 
mK/N, with A'' the total number of spins in the film) 
as a function of magnetic field H for film thickness (a) 
d = 1.0 /im and (b) d — 0.1 /im. 



FIG. 3. Transition criterion A from non-symmetric to 
symmetric phase, A (in units of mK/N), as a function 
of magnetic field H for different thicknesses d. 
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When A > 0, 5 = minimizes the energy, so 
the two minima have equal numbers of condensed 
magnons. This is the symmetric phase, with 
Nq = N^Q. When A < 0, the minimum shifts 

to = 1 — (^B-A)^ • This is the non-symmetric 
phase (S), with Nq N-q. The transition from 
symmetric to non-symmetric phase (NS) at A = 
is of the second order. There is no metastable state 
of these phases. At C = one finds S — ±Nc, 
which corresponds to a completely non-symmetric 
phase with only one condensate. The ground state 
of the non-symmetric phase is doubly-degenerate, 
corresponding to the two possible signs for S. FigjS] 
shows that for a film thickness of about 0.05 /im, 
the symmetric phase is energy favorable up to 
H = 1.2 T. For d = 0.08 fim, on increasing H to 
about 0.6 kOe, there is a transition from symmet- 
ric to non-symmetric phase. For the larger thick- 
nesses d = 0.1 /im and d — I /im, the ground state 
is non-symmetric for _ff > 0.3 kOe. 

Figjl] gives the phase diagram in {d, H) space. 
It has three different regions, separated by two 
critical transition lines, dc{H) and d*{H), corre- 
sponding to A = and C — 0, respectively. For 
d — 0.13 — 0.16 /i the system possesses re-entrant 
behavior (NS $ = tt, to NS $ = 0, to NS $ = tt) 
as H increases. As shown below, measurement of 
the contrast, or modulation depth 17J, of the spa- 
tial interference pattern permits identification of 
the different condensate phases. 
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FIG. 4. The phase diagram for different values of thick- 
ness d and magnetic field H. 

Zero Sound. — In two- condensate states the 
relative phase Scj) — — <f)- is a Goldstone mode. 
Its oscillation, coupled with the oscillation of the 
number density 5n — nq — n_Q represents a new 
type of collective excitation, which we call zero 
sound (as in Landau's Fermi liquid, this mode 
is driven by the self-consistent field rather than 
collisions). Solving a properly modified Gross- 
Pitaevskii equation (see Methods), we find its spec- 
trum. In the symmetric phase its dispersion rela- 
tion is: 

[m^ TVA^ 

V m 
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The magnon effective mass is of the order of the 
electron mass. The density of condensed magnons 
ric = Nc/V is about 10^^ cm'^ and A « 10 
mK/N. The sound speed for small k in this case 
is Vqs = \/iVcA/m, which is about 100 m/s. Near 
the transition point A = 0, the velocity of this zero 
sound goes to zero. For the non-symmetric case, 
the spectrum is: 



K + Nc(B - A){k- 1) 



fc2 



(7) 



where k 



_ {B-Af 



In the experiment of Ref. J17j , 
K ~ 10'' and^B - A = 8.4 mK/N, from which 
we estimate a sound speed of 3 x 10"^ m/s. The 
dispersions of zero sound for symmetric and non- 
symmetric cases are shown in Figjs] Note that 
the range of applicability of the linear approxima- 
tion decreases significantly for small C, where one 
of the condensate densities is small and the phase 
fluctuations grow. 
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FIG. 5. Dispersion of zero sound as a function of 
wave vector in the direction of external magnetic field 
for symmetric and non-symmetric cases, respectively. 
For the non-symmetric case, we choose H = 1 kOe and 
d = 5 /im. 

Domain Wall. — Since the ground state of 
the non-symmetric phase is doubly degenerate, it 
can consist of domains with different signs of 6 
separated by domain walls. The width w of such a 
domain wall is of the order of 



ri2 



For the 



2mNa\A\ ■ 

data of experiment |17j we estimate that w ~ 10 
/zm, and a domain wall energy per unit area of 
about 10^9 J/m2. 



Discussion 

The ground state wave function '^'{z) gener- 
ally is a superposition of two condensate ampli- 



)/vV, where 



tudes 'i'(z) = 

c±Q = V^^Qe*** and V is the volume of the 
film. The spatial structure of 4* (2;) can be mea- 
sured by Brillouin Light Scattering (BLS), with 
intensity proportional to the condensate density 
= ng -I- n_Q -I- 2^nQn^Q cos{2Qz + 
In their recent experiment, Nowik-Boltyk el al 
[T7] observed the interference pattern associated 
with the ground state. They found that the con- 
trast of this periodic spatial modulation is far be- 
low 100%; of the order 3%. The present theory 
can quantitatively explain this result. In their ex- 
periment, Ref.|17j employ d = 5.1 /xm and H = 1 
kOe. Then eq.^ for A, B and C gives A = -0.168 
mK/N, B = 8.218 mK/N and C = -0.203 mK/iV, 
so A < 0. This corresponds to the non-symmetric 
phase, where for spontaneous symmetry-breaking 
with 6 = Nq — N_q > the ratio of the num- 
bers of magnons in the two condensates is ~ 



The contrast is /3 



1*1. 



-1*1^ 



Since 



C < B and N^Q < Nq, we have /3 « 2^^^^ « 

\B-A\ ■ above values of A, B and C, j3 

is of order 2.4%, in good agreement with experi- 
ment. The smallness of C (and A) relative to B 
is associated with the large parameter d/l, where 

I — \J is an intrinsic length scale of the sys- 
tem and / ^ 10^^ m. In terms of this parameter, 



M 

B " \dl 

In the experiment of [TTl the contrast /3 reaches a 
saturation value at a comparatively small pumping 
power, corresponding to the appearance of BEC. 
This agrees with our expression for /?, which de- 
pends only on film thickness d and magnetic field 
H. By varying d and H, the contrast can be 
changed. Specifically, in the symmetric phase, 
/3 = 1; in the non-symmetric phase, /3 < 1; and 
in the completely non-symmetric case with only 
one condensate {d — d*), (5 = {). Therefore, mea- 
surement of the contrast for different values of d 
and H can give complete information about the 
phase diagram of the system, for comparison with 
the present theory. 

Figj6] plots C, A and j3 as functions of the film 
thickness d at fixed magnetic field H = 1 kOe. For 



5 



small d the system is in the high-contrast sym- 
metric state. On increasing d above dc = 0.07 
Hui, the sign of A changes, corresponding to the 
transition from the symmetric to the low-contrast 
non-symmetric phase. As d further increases, to 
d* = 0.17 /Ltm, C changes sign, and the total phase 
$ changes from tt to 0. Only at this point d* does 
the zero-contrast phase (with only one condensate) 
appear. Correspondingly, a characteristic cusp in 
the contrast /3 appears near d* . 
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FIG. 6. (a) Phase transition criterion A and inter- 
action coefficient C as functions of tliickness d for 
fixed magnetic field H — 1 kOe. (b) The contrast 
B = max imin function of thickness d for 

H = 1 kOe. S and NS denote symmetric and non- 
symmetric phase, respectively. 



To conclude, we have calculated the 4-th order 
magnon-magnon interactions in the condensate of 
a film of YIG, including magnon- non-conserving 
terms that are responsible for the coherence 
of two condensates. sFor ufhciently thin YIG 
films (d < 0.1 n) we predict a phase transition 
from symmetric to non-symmetric phase when 
the magnetic field exceeds the modest value 
of 0.2 kOe. We also predict that within the 
non-symmetric phase there is a thickness d* (H) 
where the modulation in the observed interference 
pattern should totally disappear. 

Methods 

Magnon Spectrum. — In a YIG film with an 
in-plane external magnetic field H, the magnon 
dispersion has been studied extensively . At 

low energies, YIG can be described as a Heisenberg 



ferromagnet with large effective-spin S = 14.3 [121 
24 on a cubic lattice. The Hamiltonian consists of 
three parts: 

H = -JY,^^-^j+Hd-iHY,S!, (8) 

the nearest neighbor exchange energy, the dipolar 
interaction and the Zeeman energy. We take y to 
be perpendicular to the film and the magnetic field 
to be in the plane along z. It is convenient to char- 
acterize the exchange interaction by the constant 
D = 2JSa'^ = 0.24 eVA^. The dipolar interaction 
can be calculated using the method indicated in 
Refs. [ini 130] • The competition between the dipo- 
lar interaction and exchange interaction leads to 
a magnon spectrum ujk with minima located at 
the two points in 2D wave-vector space given by 
k = (0, ±Q) (i.e. along z), with an energy gap Aq. 
For film thickness d = 5 /xm and magnetic field 
H = 1 kOe, we find that Q = 7.5 x 10"^ cm"! and 
Aq — 2.7 GHz. In the experiment of [T7] Q was 
found to be about 3.5 x 10^ cm~^, i.e. about half 
the predicted value. The reason for this may be 
associated with a rather shallow energy minimum 
as a function of wavevector. In such a situation 
small corrections to our approximate formula can 
have a large effect on the value of Q. The lowest 
band of the magnon spectrum can be calculated 
using the Holstein-Primakoff transformation [32], 
which expresses the spin operator S in terms of 
boson operators a and . 

To second order in a and a\ the Hamiltonian 
eq.Q is: 



Ho 



with 



E 

k 



(;9) 



Ak = J Ho + Dk^ + 727rM(l - Fk) sin^ + -f2TTMFk 
Bk = 727rA/(l - Fk) sin^ 6 - -/2nMFk (10) 

where Fk = (1 — e^'"^)/kd and M is the magneti- 
zation of the material {AttM — 1.76 kG). Here, 9 
is the angle between the 2D wave vector k and 
the magnetic field direction (z). Ho of eq.([9| 
is diagonalized by the Bogoliubov transformation 
Ofc = UkCk + VkcL,^ with Uk = (=%i^)^/^ and 

Vk = sgn {Bk ) ( "^2^^"'° ' -^'^^'^i'^S to the magnon 
spectrum: 



hujk 



(A 



\B, 



|2U/2 



(11) 
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Fig[T] gives the magnon spectrum along kz for typ- 
ical values of thickness d and magnetic field H. 

Number of condensed magnons Nc — 
Nq + -/V_Q. — Experimentally, the spin lattice 
relaxation time is of order 1 /is, whereas the 
magnon-magnon thermalization time is of order 
100 ns; the magnons are long-lived enough to equi- 
librate before decaying, thus making BEC possible 
After the thermalization time the pumped 
magnons go to a quasi-equilibrium state with a 
non-zero chemical potential fi. The number of 
pumped magnons Np = N{T,fi) — N{T,0), where 
NiT, Ai) = Ek e(^A.-M)/T_i : is determined by the 
pumping power and the magnon lifetime, /i is a 
monotonically increasing function of Np but can- 
not exceed the energy gap Ag. Therefore, on fur- 
ther pumping /i — Aq and some of the pumped 
magnons fall into the condensate. The equation 
Npc = N{T, Ao) - N{T, 0) thus defines the crit- 
ical line for condensation. Since Aq ^ T and 
Np <C N{T, 0) this equation can be satisfied at 
a rather high temperature. The total number of 
condensed particles is [TTJ 

TV, = Np~N{T, Ao)+N{T, 0) = N{T, ^l)~N{T, Aq). 

(12) 

In exactly 2D systems BEC formally does not ex- 
ist since in the continuum approximation the sum 
in N{T,fi) diverges. However, for strong enough 
pumping the chemical potential approaches ex- 
ponentially close to the energy gap: Aq — fi ~ 
Aoexp(-iVp/7Vo), where Nq = VTm/h?. For 
Np/No > ln(T/Ao) all pumped magnons occupy 
only one or two states ±Q. 



Eq.( 12 1 determines only the total number of par- 
ticles in the condensate. The distribution of the 
condensate particles between the two minima re- 
mains undetermined in the quadratic approxima- 
tion. To resolve this issue we have shown that it is 
necessary to consider the fourth order terms in the 
Holstein-Primakoff expansion of the exchange and 
dipolar energy. Observe that terms of third order 
occur in this expansion of the dipolar interaction, 
but since the total momentum must be zero, such 
terms vanish for the condensate momenta (0, iQ). 

Zero Sound. — We now provide details about 
calculating the zero sound spectrum. We consider 
small deviations from the static symmetric solution 
HQ — n^Q — nc/2, = TT — = 0, so that 
'^±Q — '^cl'^ + &n± with (5n_|_ = — (5n_ = (5n/2 and 



(50+ = -(50_ = (50/2. Then 
On linearizing, the energy reads: 



\ dr 



Ay 



Using the commutation relation [6(j),Sn\ = i, and 
the equation of motion ihScf) — [S(j), H], we obtain: 

H^=-J—\/Hn + AV6n, (13) 



d5n_ 



(14) 



Taking Fourier transforms of the above two equa- 
tions in coordinate and time, one arrives at the 
dispersion relations in Eq.Q. 
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